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Abstract. We present the results from a multi colour sur- 
vey performed at the Northern Ecliptic Pole (NEP). The 
survey is designed to identify the counterparts of faint 
sources from the ROSAT All Sky Survey and the IRAS 
survey and to study their optical/near- infrared proper- 
ties. We observed the central square degree around the 
NEP in the optical bands Bj and R and in the near- 
infrared band K. A shallower survey was carried out in 
the optical / band. Here we present the results of the K- 
band survey. We discuss the source counts in the mag- 
nitude range 7 mag < K < 17.5 mag and the angular 
correlation function of galaxies with K < 17.0 mag. The 
galaxy counts at the NEP display a subeuklidean slope 
in d log TV/dm. Our shallower slope does not require the 
large effects of galaxy evolution or density evolution sug- 
gested to explain the steeper slopes found in earlier sur- 
veys. The angular correlation function of galaxies follows 
a power law w\0) = A0^^ with A = (5.7±0.8) x IQ-^ and 
6 = 0.98 ± 0.15. This is in accordance with the expected 
values for stable clustering. 
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1. Introduction 

Since surveys by satellites obtained deep exposures at 
the Ecliptic Poles, the Northern Ecliptic Pole (NEP) is 
a special location in the sky. The total exposure time at 
the NEP in the ROSAT All-Sky Survey, for example, 
was < 40 000 sec as compar ed to < 400 sec at the ecliptic 
equator (Brinkmann et al. 



199£). Likewise, the deepest 



part of the IRAS survey in the far infrared (12 — 100 /im) 
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is located at the NEP, reaching flux limits down to 60 mJy 
at 60 fim (Hacking & Houck|l987|). 



We took surveys in the optical/near- infrared regime to 
identify those sources and to study their broad band 
energy distribution. Our surveys cover the central square 
degree around the NEP and were carried out in the 
optical Bj- and i?-band and in the near-infrared filter 
K. The 95% completeness limits in Bj, R, and K are 
24.0 mag, 23.0 mag, and 17.5 mag, respectively. The 
results of the identification and details on our optical 
surveys will be given in separate papers. Here we give 
the number counts of galaxies and point like sources and 
present the galaxy two-point correlation function of our 
isT-band survey, which is unprecedented as a homogeneous 
survey in dynamic range and angular size at this depth. 
Galaxy number counts and galaxy correlation func- 
tion are basic quantities in the investigation of galaxy 
evolution and galaxy distribution and have both been 
studied extensively in the optical regime. K-band selected 
samples have the advantage that the K-corrections are 
small and nearly independent of galaxy type (Poggianti 
1997 ). Moreover the absolute magnitudes of galaxies in 
K are dominated by stars of M ~ 1 Mq , which makes 
the near-infrared as a better tracer of the visible mass of 
galaxies as compared to the optical, where the absolute 
magnitudes may be strongly affected by star formation. 
With the completeness limit of K = 17.5 mag our survey 
enters a very interesting magnitude range with a high 
statistical level given by the large field. In this magnitude 
range a break in the slope of the galaxy number counts 
from 0.67 to 0.26 has been claimed by combining shallow, 
large scale surveys with deeper surveys on small fields 
(Gardner et al. 1993). From our homogeneous wide. 



medium-deep survey we can determine the exact location 
of this break. 

Our survey also bridges the gap between correlation 
functions measured in deep samples {K > 18.5 mag) 
which find amplitudes above the values expected for a 
no n-evolving t wo-point correlation function (Roche et 
al. 1998; 1999) and those valu es fro m shallower samples 
{K < 16.0 mag, Baugh et al. 1996) which are in accor- 
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Table 1. The coordinates of the NEP field centre Table 2. Completeness and classification parameters 



azooo: 18h 00m 00.0s Szooo- 66° 33' 38.6" 
hi: 96°.38 bir. 29°.81 



dance with stable clustering (no evolution). 



Bin 


Area [de 




iVtot 




K < 16.0 mag 


0.93 de^ 




3572 


4 


16.0 mag < K < 16.5 mag 


0.93 dei 




1430 


10 


16.5 mag < K < 17.0 mag 


0.91 de^ 




2080 


120 


17.0 mag < K < 17.5 mag 


0.61 de^ 




1733 


371 



2. Observation and data analysis 

2.1. Observation and photometric calibration 

The observations were made with the 2.2-m telescope on 
Calar Alto, Spain, during two observing runs from July 22 
to 26, 1994 and August 8 to 10, 1995 using the NIR cam- 



era MAGIC (Herbst et al. 1993). This camera is equipped 
with a NICMOS 3 array, similar to the system described 
by Hodapp et al. (|99^ with 256 x 256 HgCdTe diodes. 
In the Wide Field Mode, chosen for our project, the field 
size is 6.9' x 6.9' with a scale of 1.6"/pixel. 
The coordinates of the NEP in various coordinate systems 
are given in Table 0. To cover the central square degree 
around the NEP, we observed a regular grid of 9 x 9 fields. 
On each field several images were taken on 16 different po- 
sitions with non-integer pixel offsets in right ascension and 
declination between the individual positions. This dither- 
ing allows for averaging potential variations of the re- 
sponse on a sub-pixel scale as well as for corrections of in- 
dividual hot pixels, both being essential in undersampled 
imaging. In a photometric night we carried out a snapshot 
survey where every field was observed in two positions 
to obtain homogeneous photometry. Including the snap- 
shot survey the total integration time per field is 648 sec. 
The data were obtained using a K'-tilter (Wainscoat & 
Cowie 1992) with central wavelength Xcentrai = 2.10 /im 



and AA = 0.34 fj.m width (Herbst |1995|). 
In a photometric night we observed NIR standard stars 
from the lists of Elias et al. ( 1982| ), Zuckermann & Beck- 
lin (1987), and Casali & Hawarden (1992) and used the 



i^-magnitudes to determine zero point and extinction. 
Within the range —0.036 mag < H — K < 0.325 mag cov- 
ered by the standard stars, the zero point does not depend 
on colour va H — K. 

Flux calibration was established using the snapshot survey 
taken within three consecutive hours in a single photomet- 
ric night. These data were reduced and coadded indepen- 
dently from the main survey to obtain a homogeneous set 
of photometric zero-points. This was transferred to the 
deeper, coadded survey data for each field individually, 
using several stars in every case. 

2.2. Standard reduction 

Standard reduction included dark subtraction, sky sub- 
traction, and fiatfield division of the individual images. 
Dark exposures were taken during daytime with the same 



camera settings that were used for the scientific exposures, 
and flatfields were constructed from the differences be- 
tween images of an illuminated and dark screen in the 
dome. For every field a sky frame was constructed by scal- 
ing the dark corrected exposures of this field to the same 
level multiplicatively and removing the objects by com- 
puting the median for each pixel. 

The spatial offsets of the images taken on different posi- 
tions within one field were determined using bright stars. 
Before coadding the images were resampled to a grid with 
0.4" pixels. Due to the resamphng it was possible to coadd 
the images on a finer grid. This resulted in a better def- 
inition of the objects shape and a stable performance of 
the source-detection routines. 

Further analysis was then done only on the central area 
of each field which was observed with the full integra- 
tion time. This leaves an effective survey area of 0.93 deg^, 
spread over the full square degree with small gaps between 
adjacent fields. 

2.3. Object detection and photometry 

The detection of sources was carried out with FOCAS 
(Valdes |1994| , Jarvis & Tyson [1981D . We tested the refia- 
bility of the detection process as a function of the input 
parameters on three fields which were observed with an 
integration time of 2 160 sec, more than three times the 
amount spent on the other survey fields. In these fields 
we coadded the total of all images to deep frames as well 
as subsets of images to frames of integration times and 
depths equal to those of the entire survey. Then we per- 
formed source detection on all frames with various sets 
of detection parameters. The reliability of the detection 
was analyzed by identifying the objects detected on the 
shallower frames in the deeper ones. The detection pa- 
rameters for the survey were chosen such that only < 1 % 
of the sources in the shallow frames could not be found in 
the deep frames, which gives a reliability of > 99 % for the 
sample of the detected objects. This high reliability was 
achieved by accepting objects above 16 a only. 
FOCAS performs different types of photometry: the 
isophotal flux Li is measured within the detection area, 
the so called total flux Ltotai, which is determined in the 
area expanded by a factor two around the detection area 
and the flux measured within a specifled aperture around 
the object, I/;/ca- Ltotai has the advantage to give good 
flux estimates independent from the extent of the objects. 
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Fig. 1. The number density for extended and point-like 
sources, labeled as galaxies and stars, respectively. 



because the flux is determined in an area adapted to ob- 
ject sizes. However, Ltotai misses flux at the faint level 



(see e.g. Saracco et al. 1997). To take this into account we 
chose Ltotai for the bright sources and Lifca, measured in 
an aperture of 7.2" diameter, for the faint sources. The 
final value depends on the size of the object. If Ltotai is 
determined in an area smaller than 41 arcsec^ (the area 
corresponding to the diameter of the aperture), i;/ca is 
used instead. The transition from Ltdai to L;/ca occurs 
ai K ^ 16.5 mag. We compared the magnitudes of the 
same objects found on the set of shallow and deep frames 
mentioned above to confirm that the magnitude differ- 
ences are caused by statistically distributed measurement 
errors and did not detect any systematic trend. 

2.4- Morphological classification 

Since the point spread function (PSF) of the J^-images 
is grossly undersampled, we made no attempt to clas- 
sify the objects on the basis of the i^T-data. Instead, 
the objects were classified on the basis of the morpho- 
logical identifying of their counterparts in the i?-survey 
which has an effective PSF of ^ 1.2". In the i?-survey (of 
which details are given in Kiimmel 1998 and Kummel & 
Wagner in prep.) the classification is based on the FO- 
CAS resolution classifier (Valdes 19821 ). The classification 
into point-like and extended objects is reliable for objects 
R < 21.5 mag. Towards fainter magnitudes extended ob- 
jects arc progressively misclassifled as point-like objects 
because their extended parts vanish in the background 
noise. In the R band we use a statistical classification 
by extrapolating the number-counts of point-like objects 
(dlogiV/dm = const) into the range of unreliable clas- 
sification to estimate the fraction of extended objects. 
The (small) fraction of sources detected in K whose R- 
band counterparts are too faint for a reliable morphologi- 
cal classification were divided into point-like and extended 
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Fig. 2. The galaxy counts from the NEP compared to the 
results of other surveys covering areas on the square degree 
scale. The solid line shows predicted counts from Jimenez 
& Kashlinsky ( |1999D (their flo 0.2, = 0.0 model). 



sources based on the fraction of these types at the corre- 
sponding magnitudes of the i?-band counterparts. This is 
strictly valid only if point-like and extended objects have 
similar R-K colour for R < 21.5 mag 

2.5. Completeness 

The completeness was determined using the set of 
deeper images with longer integration times, mentioned 
in Sect. p.3| which were used already to find the optimum 
set of detection parameters. The object lists of those im- 
ages were divided into those found in the shallower images 
and those found only in the deep images. This ratio was 
fitted by a completeness function f{mag) 



f{mag) = 



exp 



/ mag — mag^Q 
[ b 



1 



(1) 



with mag^Q being the magnitude where 50 % of the ob- 
jects detected in the deeper images have been found in 
the shallower images. The parameter & is a measure of the 
width Amag = mag2—magi where completion drops from 
f{magi) = 1 — e to f(rnag2) = e. We fitted f{mag) to a se- 
ries of deep/shallow matches taken under different condi- 
tions with different image quality. While h was found to be 
independent of image quality {h = 0.26), magso varies by 

0. 3 mag, reflecting the basic parameters of image quality, 

1. e. background noise and FWHM of the PSF. Establish- 
ing the relations between those parameters and mag^o and 
using these relations and the corresponding zero points we 
computed mag^o for every individual survey image. We fi- 
nally set a completeness limit magcompi = ™a550 ~ 0.5. 
Down to this level f{mag) was found to be > 0.9 (i.e. 
indistinguishable from 1.0) in our deep fields. We do not 
use f(mag) to correct counts fainter than magcompi in this 
study. 
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Table 3. Quantitative comparison of galaxy counts 



Survey 


area 


slope 


amplitude 


range 


NEP„ 


0.93 


0.56 ±0.01 


97 ±4 


10.5-17.0 


NEPft 


0.93 


0.41 


175 


16.5-17.5 


Szokoly 


0.6 


0.56 ±0.02 


98 ±5 


12.75-16.75 


Huang 


9.8 


0.65 ±0.01 


103 ±2 


11.0-16.0 


Gardner 


9.8 


0.63 ±0.01 


103 ±2 


10.0-16.0 



3. Source counts 

3.1. Star and galaxy counts 

Fig. |l| displays the object counts for point-like and ex- 
tended source, labeled as stars and galaxies, respectively. 
No attempt was made to correct object counts beyond 



magcompi as defined in Sect. 2.5. The counts were derived 
only from fields complete to the specific magnitude. The 
area of the sub-surveys complete to various depths are 
given in Table ^ Furthermore, Table |^ shows the number 
of objects with statistical classification in R, Nstat, and 
the total object number in that bin, iVtot. 
The star counts are in good agreement with the theoretical 
expected values from Wainscoat et al. (1992), calculated 
for III = 90°. 0, bji = 30°. 0, which is in the vicinity of the 
NEP (see Table 

For K < 17.0 mag the galaxy counts are well approxi- 
mated by a power law with d log N/dm = 0.56. The counts 
in the faintest bin differ significantly from this slope and 
the two faintest bins result in d\ogN/dm = 0.41, estab- 
lishing a break or the beginning of a roll-over in the range 
16.5 < K < 17.0. Both the bright and the faint end slope 
are below the Euclidean value of dlogiV/dm = 0.6. 

3. 2. Interpretation of the galaxy number counts 

To compare our galaxy counts with the numbers from 
other surveys covering areas > 0.6 deg^. Fig. |^ gives our 
data together with the counts from Gardner et al. ( 1996D , 
Huang et al. ( |1997D and Szokoly et al. ( |1998D . In Fig. | 
the Euclidean slope 0.6 has been subtracted from the log- 
arithm of the counts in order to expand the ordinate and 
to make differences between the counts clearly visible. 
We fitted power laws of the form 



N{mag) = a * IQbH^nag-lS.Q) 



(2) 



to the counts plotted in Fig. || and give the parameters 
of those fits as well as quality estimators in Table ^. The 
first two columns give the reference and the survey area 
(in deg^). Then the slope b, the amplitude a, and magni- 
tude range of the fits are listed. The counts at the NEP 
are split into two regimes, referenced as NEPa and NEPf,, 
because counts at the faint end (NEP;,) do not follow the 
extrapolation of the power-law of the K < 17.0 mag ob- 
jects (NEPa in Table I). 

The normalization of the power-laws for the different 



galaxy counts with K < 16.5 mag, listed in Table ^, agree 
well among each other within the errors. The slopes of 
the power laws differ. While the surveys of Gardner et 
al. (1996) and Huang et al. (1997) find values around 
b sa 0.64, the galaxy counts of Szokoly et al. (1998) 
and at the NEP both result in the subeuklidean value 
b — 0.56. The difference between those two groups can 



be seen clearly at the counts in Fig. 
the data points from Gardner et al. 



Power-law fits to 



996) or Huang et 



al. ( 1997 ) have a positive slope, which means a supereuk- 
lidean slope for the p ure c ounts, whereas fits to the data 
from Szokoly et al. ( 1998 ) and the NEP have negative 
slopes corresponding to a subeuklidean slope for the pure 
counts in dlo gN/d m. 

Huang et al. ( 1997 ) explain the large slope of their counts 
with a local under-density of galaxies. Quantitatively they 
determine that the normalization of the luminosity func- 
t ion d oubles its value from z = to 0.2. Phillips & Turner 
(1998) analyze all fsT-band counts to K = 17.0 mag known 
to mid- 1997 and find two possible explanations for the 
steep slope: Either an evolution of the normalization of 
the galaxy luminosity function by a factor 1.7 — 2.4 to 
a redshift of 0.1 — 0.23 or an unexpectedly strong low 
redshift evolution in K which leads to evolutionary and 
cosmological corrections as much as 60 % larger than ac- 
cepted values. 

Our shallow, subeuklidean slope of the galaxy counts at 
the NEP argues against a local underdensity of galaxies or 
strong luminosity evolution. The value dlog A^/drn = 0.56 
is even smaller than the range 0.602 — 0.623, which were 
computed in Huang et al. (1997) for different luminosity 
functions and cosmological corrections based on constant 
galaxy luminos ities. 



Gardner et al. ( 1993 ) reported a break in the slope of the 
galaxy number counts from dlogiV/dm — 0.67 to 0.26. 
That break was claimed to occur between K « 15.0 mag 
and K « 18.0 mag. Our data at the NEP show with high 
statistical significance that this break in slope is less ex- 
treme, and starts at K 16.5 mag with a flattening to 
0.41 at K Ri 17.5 mag. A further flattening to slopes < 0.3 
must occur at magnitudes K > 17.5 — 18.5 mag to match 
number densities to deeper source counts from e.g. Mous- 
takas et al. ( 1997 ) and Djorgovski et al. ( 1995 ). 
Subeuklidean counts in the magnitude range K — 15.0 — 
17.0 mag are supported by the recent compilation of 
galaxy evolution models from Jimenez & Kashlinsky 
( 1999| ) . In their model the density of galaxies is assumed to 
be independent of redshift. Their predicted i^-band counts 
(for Hq = 0.2 and SIa — 0.0) are shown with the solid line 
in Fig. 1^. Although there is an offset in the normalization 
of ~ 0.15 dex with respect to all data displayed in Fig. ^ 
the slope of the predicted counts drops below euklidean 
values at K > 14.8 mag. 



4. Angular correlation function 
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Table 4. The angular correlation function for the various samples 
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Sample 


selection 


Number of 


A free, 5 fixed 


A and S free 




K [mag] 


R-K [mag] 


galaxies 


A [10-^] 


S 


a^/A 


A [10-^] 


S 


a^/A 


12.0 - 16.0 


all 


491 


9.97 ± 2.77 


0.8 


1.83 


2.78 ± 3.42 


1.14 ± 0.27 


2.84 


16.0 - 17.0 


all 


1234 


5.75 ± 1.07 


0.8 


1.84 


2.33 ± 2.13 


1.02 ± 0.21 


2.44 


12.0 - 17.0 


all 


1725 


5.70 ± 0.77 


0.8 


1.83 


2.77 ± 1.78 


0.98 ± 0.15 


2.30 


12.0 - 17.0 


>3.49 


858 


9.23 ± 1.55 


0.8 


1.85 


3.12 ± 2.47 


1.08 ± 0.18 


2.81 


12.0 - 17.0 


<3.49 


867 


7.30 ± 1.52 


0.8 


1.83 


3.48 ± 3.50 


0.99 ± 0.23 


2.32 



4-.1. The calculation of the angular correlation function 
The angular correlation function w{9) is defined by 

dP ^ N^[l + w{e)]dQidn2, (3) 

where N is the mean density of galaxies per steradian and 
dP the joint probability of finding one galaxy within the 
solid angle df22 in the projected distance 9 away from 
a galaxy within the solid angle dfii (Groth & Peebles 



1977). w{9) is determined by assembling the normalized 
histograms of logarithmic distances for the data sample, 
{DD), a random sample, (RR), and for the d ata-ra ndom 
sample, (DR), and applying (Landy & Szalay 1993 ): 



w{9) = 



{DP) - 2{DR) + (RR) 
(RR) 



(4) 



The measured function w''^{9) computed in Eq. (|^) is bi- 
ased towards low values with respect to the true function, 
w\9), by 



w"'{9) =w\9) 



(5) 



(Baiighet al. |l996i Roche & Bales |1998| , Lidman & Peter 



1996 ), but this bias ct^ can be estimated as 
^ / w\9i2)d^idn2 



(6) 



(Groth & Peebles 1977). The angular correlation function 
is usually described by a power-law w^{9) = A9^^ . If the 
exponent 5 is known, the amplitude of w*{9) is determined 
by fitting the measured correlation to the function 



w"'{9) =A{9-'' - a' I A). 

If the value for b is not known, Eq. 
to be iterated. 



(7) 

and Eq. (0) have 



^.2. w{9) of galaxies K < 17 mag 

We calculated the angular correlation function of galaxies 
brighter than K = 17 mag. To avoid spurious detection of 
the filaments around the planetary nebulae NGC 6543 a 
sufficiently large region was masked out, leaving a sample 



of 1 725 o bjects in an area of 0.9 deg . Following Baugh et 
al. (1996) we use 



5w"(6i) = 2 * y/{l+w"'{9))/[DD]) 



(8) 



to derive errors, were [DD] denotes the non normalized 
histogram of (DD). 

Table ^ give the results of fitting w™ (9) to Eq. (|^) in case 
for a fixed d = 0.8 and for 6 being a free parameter. This 
was computed for the total sample as well as for subsam- 
ples split at K = 16 mag. 

Fig. ^ shows the bias-corrected angular correlation func- 
tion of the total sample with the power law fitted to the 
data {5 = 0.8) on scales 10" < 6* < 33'. While the lower 
limit is dictated by the small number of pairs in [DD], 
the upper limit is set to half the largest size of the sur- 
vey field. The measured correlation follows the power law 
very well on all scales. There is no indication of a break 
or discontinuity on the scale of the distance between two 
adjacent fields (marked by the arrow in Fig. |3|). Even the 
largest deviation at the widest angular scale {9 = 0.44°) 
is at the level of 2.5 a, i.e. not significant. 
With a power-law slope oi S ^ 1.0, and range of more 
than two orders of magnitude in 9, a strong variation of 
the surface density about the mean value of 0.5arcmin~^ 
is implied. Due to the still low average surface density, 
the peaks in the surface density distribution cannot be at- 
tributed to clusters in an unambiguous way. 
We investigated whether the angular correlation depends 
on galaxy colour. The total set of objects is divided into 
a red sample with R ~ K > 3.49 mag and a blue sample 
with R — K < 3.49 mag. The two subsamples are selected 
to include approximately the same number of objects (see 
Table ^). The angular correlation function was computed 
in the same way as for the total sample and values for the 
amplitude and scale-length derived from the fits to the 
correlation functions are given in Table ^. 
Fig. 1^ shows the bias corrected correlation functions of the 
blue and the red sample in comparison to the total sample. 
As can be seen in Fig. || and in Table ^ within the errors 
the blue galaxies show the same clustering properties as 
the red galaxies. 
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Fig. 3. w*{9) for the K < 17.0 mag and for blue and the 
red subsamples. Abscissa values of the red and the blue 
sample were offset by Alog(6') = ±0.02. The solid line 
displays the power law fitted to the total sample {S = 
0.8) and the arrow indicates the separation between two 
adjacent fields 

4-3. The angular correlation function of K -selected 
galaxies 

The two-point correlation function is approximated by a 
power-law with two free variables. It is important to note 
that these are not derived independently. The comparison 
of the amplitude of the correlation is subject to differ- 
ences in the angular range involved in the fitting even for 
fits with fixed slopes S. 

Fig. ^ gives the change of amplitudes of the angular corre- 
lation function with limiting magnitude in K. The values 
determined in Sect. ^ are displayed together wi th all other 
measurem ents k nown to us (f rom Baugh et al. 1996 , Carl- 
berg e^l. |l997t Roche et al. |l998i |l999| ). From Baugh et 
al. ( 1996 ) two values for the amplitude are given in Fig. 
^ since the amplitude obtained in two different fields are 
discrepant at the bright end. 

The lines in Fig. ^ show the theoretically expected values 
for the amplitude of the angular correlation function as 
computed in Roche et al. ( 199^ 199E ) for various mod- 
els of the spatial two-point correlation function. Details of 
those models can be found in these papers and references 
therein. 

The amplitude of our K < 16.0 mag galaxies is in accor- 
dance with the values from Baugh et al. (1996) and does 



not exclude any of the models A to E due to the large 
error. The amplitude at the faint level is in good agree- 
ment with the values expected from models A and C. The 
difference between the data and models B and D is 2.2 a, 
thus these models cannot be rejected. Model E is rejected 
at the 3.3 a level. We want to point out that the ampli- 
tudes compared in Fig. 5 are derived assuming S = 0.8. 
Piecewise fits of w"^{9), shown in Fig. 3 reveal that nar- 
row angle surveys would obtain lower normalization and 



Fig. 4. The w{9) amplitude of galaxies in our survey at 
the NEP, compared with amplitudes from the surveys of 
Baugh et al ( |1996| ), Roche et al. ( |1994 |1999D and Carlberg 
et al. (1997). The lines show the theoretical models from 



Roche et al. (199?, 199S). 



wide angle surveys larger ones. Comparison of dA/dK as 
shown in Fig. 5 rely on w"^{9) being well described by a 
single power-law. 

Our data favour model A and C, which were computed as- 
suming stable clustering and a significant higher normal- 
ization of the clustering strength for early type galaxies, 
as measured by Guzzo et al. (1997). 

A stronger clustering of early type galaxies should result 
in a larger angular correlation amplitude for red galax- 
ies. This is in contradiction to the results from Sect. 4.2 
where the amplitude of the red sample (early type galax- 
ies) , was not significantly larger than the amplitude of the 
blue sample (late type galaxies). This disagreement could 
be caused by the diluting effects of the if -band selection 
of the sample which results in a bro ader distribution of the 
red galaxies in dN / dz (Roche et al. 1999) and a decrease in 
the angular correlation amplitude of the red sample. This 
would still imply the existence of several nearby clusters 
with increased numbers of early-type galaxies, which are 
not detected. 

Finally, we want to stress that further confirmation for the 
high correlations (see Fig. ^ found in K > 19.0 mag sam- 
ples is needed despite existing measurements of the angu- 
lar correlation function in samples down to K = 21.5 mag 



(Carlberg et al. 1997). Due to the small sample sizes used 
in these studies, the derivation of amplitudes (four times 
smaller than those found at K<17 in the present study) 
are dominated by Poisson errors and uncertainties intro- 
duced by the correction for stellar contamination. The 
measurement of the small signals expected at K 19 mag 
would require samples of several 10** sources which are not 
yet available. 
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